AFCRL-65-440 
JUNE  1965 

PHYSICAL  SCIENCES  RESEARCH  PAPERS,  NO.  114 


O 

>s 

\ 

a 

<1. 


SPACE  PHYSICS  LABORATORY 


PROJECT  8606 


AIR  FORCE  CAMBRIDGE  RESEARCH  LABORATORIES 


L.  G.  HANSCOM  FIELD,  BEDFORD,  MASSACHUSETTS 


Distor^on  of  a  Magnetic  Field  by  the 
Motion  of  a  Cylindrical  Conductor 

JOSEPH  G.  KEUEY 
RONALD  S.  FINN,  .1/LT,  USAF 


OFFICE  OF  AEROSPACE  RESEARCH 
United  States  Air  Force 


Abstract 


Motion  of  a  conductor  relative  to  a  magnetic  field  distorts  the  field.  This 
paper  considers  a  cylindrical  slug  moving  in  a  two-dimensional  magnetic  field, 
represented  by  the  vector  potential  Aq  9,  where  3  Aq/9  0  =  0.  Maxwell's  equations 
are  solved  for  the  distorted  potential  A  6  in  the  form  of  a  rapidly  converging 
series  A  =  Z)  A^^.  The  A^'s  are  given  in  a  form  suitable  for  evaluation  by  a  digital 
computer.  The  nonequivalence  of  the  apparently  analogous  problem  of  a  stationary 
slug  in  a  time-varying  field  is  noted  and  discussed. 
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Distortion 


of  a  Magnetic  Field  By  the  Motion 
of  a  Cylindrical  Conductor 


1.  INTHOniCTION 


1.1  Itac'kpround 

By  Faraday's  law  a  change  of  the  nr.agnetic  field  in  a  conducting  slug  induces 
currents  in  that  slug.  This  change  may  be  either  time  variixtion  of  the  field  on  a 
stationary  slug  or  motion  of  the  slug  in  a  constant  field.  By  Ampere's  law  these 
currents  produce  an  induced  magnetic  field  whicn  perturbs  the  original  field.  Such 
distortions  have  been  use...  to  measure  the  conductivity  of  various  moving  (Lin, 
Resler  and  Kantrowiiz,  1955)  and  stationary  (Chambers  and  Park,  1961)  slugs. 
Theoretical  analysis  has  been  limited  to  a  rectangular  slug  in  a  one -dimensional 
field  (Oddson,  1963).  These  devices  have  been  calibrated  by  measuring  slugs  of 
known  conductivity.  However,  some  phenomena,  such  as  skin  effect,  which  were 
not  a  problem  in  these  cases,  can  be  difficult  to  duplicate  in  the  calibrating  slugs. 
To  extrapolate  the  calibration  to  these  cases,  one  must  solve  the  problem  for  a 
cylindrical  slug  in  a  two-dimensional  field. 

1.2  Prohlcm 

Consider  an  initial  field  represented  by  the  vector  potential  Aq  ^  with  the 
condition  that  3  Aq/3  0  =  0.  A  cylindrical  slug  of  radius  b,  length  2jE,  and 


(Received  for  publication  6  April  1965) 


2 

conductivity  o,  moves  into  this  field  at  velocity  v.  (Figure  1.)  The  slug's  midpoint 
is  at  position  z  =  vt  at  time  t.  What  perturbation  in  A  will  be  caused  as  a  result  of 
this  motion?  Since  this  problem  involves  spatial  variations  of  the  field,  it  will  be 
referred  to  as  the  Space  Case  or  Space.  Consider  also  the  problem  of  an  identical 
but  stationary  slug  with  its  midpoint  at  the  point  Zq  -  const,  in  an  initial  field 
represented  by  the  vector  potential  Aq  (t)  Again  with  the  condition  that 
9  Aq/9  6=0.  What  perturbation  in  A  w'ill  be  caused  as  a  result  of  this  time  varia¬ 
tion  of  the  initial  field?  Since  this  problem  involves  temporal  variations  of  the 
field,  it  will  be  referred  to  as  the  Time  Case  or  Time. 


Figure  1.  Geometry  of  Problem.  Acylindricalslugof  Length  2f  and  radius 
b  moves  along  the  z  axis  with  velocity  v.  At  time  t  the  midpoint  of  the  slug  is  at 
the  point  where  z  =  vt.  A  current  f^  flows  in  a  single  turn  loopat  r^,  Zq=  0 


For  the  purpose  of  this  discussion  let  Aq  be  created  by  current  Iq  flowing  in 
the  loop  at  r^,  Zq  -  0.  More  complex  fields  of  this  symmetry  may  be  represented 
by  a  number  of  such  loops. 


3 


i 


1.3  Method  of  Solution 

The  symmetries  postulated  above  allow  one  to  reduce  Maxw?ll's  equations  for 
moving  media  to  one  scalar  differential  equation  for  A,  ^vhere  A  is  the  vector 
potential  for  the  total  magnetic  field,  both  initial  and  induced.  Let  A  equal  Z)  A^^ 
vdiere  A^  is  the  order  perturbation.  Then,  by  physical  reasoning,  Aq  generates 
first  order  eddy  currents  which  create  the  first  order  perturbing  field  A^. 
Similarly  potential  A^  generates  currents  which  create  potential  Ag,  and  so  on. 
This  can  be  carried  to  as  many  orders  as  accuracy  requires.  Rationalized  mks 
units  will  be  used  initially.  Later  all  quantities  will  be  made  dimensionless. 

2.  SOLUTION 

2.1  Malhcmatical  Lxprcssion  of  Problem 

Consider  Maxw-ell's  equations  for  moving  media  (Panofsky  and  Phillips,  1955). 

V  .  E  =  0 

V  •  B  =  0 

V  =  /‘oio 

V  X  B  =  Pq  Jq  +  fiQOE  +  (0  E/dt) 

V  X  E  =  -  a  B/a  t 

Jq  is  the  current  density  which  produces  the  initial  field. 
o  is  the  conductivity  of  the  medium, 
y  is  the  velocity  of  the  medium  relative  to  the  fields. 

Charges,  both  true  and  those  due  to  polarization,  are  assumed  to  be  nonexistent. 
The  velocity  of  the  medium  is  much  less  than  the  velocity  of  light.  Maxwell's 
equations  outside  of  the  medium  are  obtained,  if  the  velocity  and  conductivity  are 
set  equal  to  zero. 

The  fields  may  be  expressed  by  the  vector  potential  A. 

E  =  -a  A/at 
B  =  VX  A  . 


(Space) 

(Time) 

(1) 


(2) 
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Combining  Eq.  (2)  and  Eq.  (1),  one  finds  that 

VX  (VX  A)  =  Hq  Jq  -  (3  A/3t)  +  Mq  ^  v  X  (VX  A)-  Mq^q  (Space) 


VX  (VX  A)  =  MoJo  ■  (9  A/3t)  -  fi^e^  (3^A/3t^).  (Time)  (3) 

F._r  this  problem,  the  vectors  have  the  following  forms: 

io  "  ‘io  ^ 

V  =  V  k  (4b) 

B  =  r  +  k  (4c) 

E  =  E  ^  (4d) 


due  to  the  symmetries  postulated  above. 

From  Eqs.  (2),  (4c)  and  (4d),  it  follows  that 

A  =  A  0  (5) 

since  it  has  been  assumed  that 

3  E/3  0  =  a  B/3  0  =  0  (6) 

then 

3  A/3  0  =  0.  (7) 

From  these  considerations  it  can  be  shown  that  only  the  0  component  of  Eq.  (3)  is 
important. 

(3/3r)[r  Ma/c^r)  (r  A)]  +  (3^A/3z^)  =  -  o  v  (3  A/a z) 

+  Mq  (3A/3t)  +  Uq  (3^A/3t2)  (Space) 

(3/3r)  [  r  ^  (3/3r)  (r  A)  j  +  (3^A/3z^)  =  -fx^  +  Mg  (3A/3t) 

+  Mgeg  O^A/at^) 


(Time)  (8) 


i 
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-2  2  2 
Note  that  =  c  ,  where  c  is  the  velocity  of  light.  The  term  '9  A/9z  )  is 

of  the  same  order  as  v  ^  (9^A/9t^)  (Space)  or  0)^  ^  r  ^  (9^A/9t^)  (Time).  The 

frequency  a)«  is  some  fundamental  frequency  of  the  Time  Case  so  chosen  that 

^  2  2 
Wq  Cq  (Time)  is  on  the  order  of  v  (Space).  Thus  the  term  Mq  ^0  A/9t  )  - 

(v/c)^  V  ^  (9^A/9t^)  (Space)  =  (cOq  ^q/c)^  cUq  ^  r^  ^  (9^A/9t^)  (Time)  may  be 
omitted. 

This  makes  Eq.  (8)  a  diffusion  equation.  Since  physically  this  is  a  diffusion 
problem,  the  approximation  is  valid. 

All  quantities  will  now  be  made  dimensionless  as  follows:  Choose  a  charac¬ 
teristic  current  and  a  characteristic  length  r„.  Multiply  each  term  in  Eq.  (8) 

2  _i  u  u 

by  r^  (/1q  Iq)  .  Substitute  dimensionless  terms  for  the  resulting  ratios  as  indi¬ 
cated  in  Table  1.  For  the  remainder  of  this  paper  only  these  dimensionless  terms 
will  be  used.  The  dimensionless  equivalent  of  Eq.  (8)  is 

(9/9  r)  [r  ^(9/9r)(rA)]  +  (9^A/9z^) 

■  "Jq  +  Q!  [  (9  A/9z)  +  (9A/9t)]  (Space) 

=  -Jq  +  01  [  9A/9t]  .  (Time)  (9) 


Table  1.  Relationship  BetA'een  the  Ratios  of  Quantities  in  mks  Units  and  the 
New  Dimensionless  Quantities 


Ratios 

Dimensionless  Quantities 

r/ro 

r 

z/ro 

z 

vt/r^  (Space) 
oJq  t  (Time) 

t 

> 

o 

o 

A 

2/, 

>io  ''o  ^^0 

jo 

o 

u 

b 

S.lv 
'  *^0 

£ 

o 

u 

o 

u 

1 
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Table  1.  Relationship  Between  the  Ratios  of  Quantities  in  mks  Units  and  the 
New  Dimensionless  Quantities  (Cont) 


Ratios 

Dimensionless  Quantities 

^o/^o 

==0 

Hq(T  V  Vq  (Space) 
Mq  ^  *^0^  (Time) 


2.2  Functions  to  he  L  scd 

Consider  first  the  potential  a  due  to  a  current  I  in  a  loop  r  ,  z  .  The 
^  n  n  n  n 

cu'"rent  density  is  given  by 

6  (z  -  z^)  6  (r  -  r^).  (10) 


The  Dirac  delta,  6  (y  -  yj^)/  is  defined  by  Eq.  (11). 


yi 


5  (y  -  y^)  d  y  =  1  yi  <  y^  <  yg 

=  0  otherwise 

g  (y)  5  (y  -  y^)  d  y  =  g  (y^)  y^  <  y^  <y2 
=  0  otherwise. 


The  potential  is  given  by  (Stratton,  1941) 

a„  =  (2«)‘"  I  (r  /r)^/2  V 

n  n  n'  n 

'’’he  following  definitions  are  required: 

^n  =  /^^^[l  -  X  sin^9l^/^  d  (> 

0 


(11) 


(12) 


(13) 
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and  are  the  complete  elliptic  integrals  of  the  first  and  second  kinds  (Jahnke 
and  Emde,  1945).  The  function  f  (x^^)  is  shown  in  Figure  2;  a^  (r^,  z^)  is  mapped 
in  Figure  3. 

Azimuthal  currents  in  a  solid  may  be  regarded  as  a  continuous  distribution 
of  such  rings.  To  obtain  the  potential  due  to  such  currents,  one  must  integrate 
over  all  such  rings  and  divide  by  an  appropriate  normalizing  factor.  Thus, 


b  t+£ 

A  (r,  z)  =  fn  f.f,a  dz  dr 
n  JQ  J  t-£  n  n  n 


z  b  £ 


L  r  a  dz  dr 
•'0  J  A  n  n  n 

^0  ^ _ 

z  b  £ 


(Space) 


(Time)  (14) 


It  should  be  noted  that 


^n  =  ^n  ^n>- 


(15) 


Since  the  normalizing  factor  will  later  cancel,  large  values  of  £  present  no  problem. 
Note  that  n  refers  to  the  n'”-order  perturbation.  The  rings  are  not  discrete  and 
therefore  are  not  numbered 

One  more  function  is  needed.  We  define  A  ,  to  be  the  function  A  ,  with 

n,  n-1  n-l 

r„  z  substituted  for  r,  z.  Similar  definitions  may  be  formed  for  x  , ,  K  , , 
n,  n  n,  n-l'  n,  n-l' 

and  a^  is  the  potential  at  r^,  z^  due  to  a  current  in  the  loop  at 

n-l  n-l 

Since  a^  is  the  potential  of  a  current  loop,  it  must  satisfy  an  equation  similar 
to  Eq.  (9)  except  a  =  0. 


(a/9r)  [r  ^  (9/Or)  (ra^,)]  +  (9^  a^/9z^)  =  -1^6  (r  -  6  (z  -  z^).  (16) 

Since  a  differential  with  respect  to  r  is  unaffected  by  integration  over  r  , 

n 

Eqs.  (14)  and  (16)  show  that 
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(a/8r)[r'^  (0/9r)  (r  A^)]  +  (3^  A  fdT?) 


.b  +.+£ 


=  -{2btl  ,r  .  r^,5  .  7.„)  dr„ 


(Space) 


-  -tZb£)  ^  f  /  °  *  I  (r  z  )6  (r  -  r  ) 
0  Zq-£  n  n  n  n 


6  (2  -  Zjj)  dz^  dr^{Time)  (17) 


■Jsing  the  definition  of  6  (y  -  y  )  in  Eq.  (11), 


one  may  show  that 


(9/3r)  [r  ^  (3/ar)  (r  A^)]  +  (3^  Ajdz?) 

=  -(2  b£)  ^  (r,  z)  r  <  b,  t  -  £  <  2  <  t  +  C 


"O  '  ^  1  ^'  <  2o  +  f 


otherwise . 


2.3  Soluiion 


Let  the  solution  to  Eq.  (9)  be 
A  =Ea^. 


Tbe  A^'s  will  now  be  evaluated. 
Outside  of  the  slug,  Eq.  (9)  becomes 


(9/9 r)  [r"  (a/9r)  (r  A)]  +  (9^A/dz^)  =  . 


(Space) 

(Time) 


If  Eq.  (19)  is  used  to  substitute  for  A  on  the  left-hand  side  of  Eq.  (20),  the  A 
term  will  contribute  -Jq  to  the  right-hand  side  and  the  other  .A  >s  will  contribute^ 
nothing  to  the  right-hand  side  and  the  solution  is  obviously  valid.  Inside  the  slug 
Eq.  (9)  can  be  satisfied  iij  the  same  manner  if 


Of  [(9A/9z)+  OA/at)]  =  -t2bf)‘'^  2  I  (r  z) 

rv  '  '  • 


(21) 
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If  Eq.  (19)  is  now  used  to  replace  A  on  the  left-hand  side  of  Eq.  (21) 


(-2b£)‘^  E  I  (r.  z)  = 
n=l  ” 


CO 

0!  E  [  (y  A^_j^/3z)+(3  A.j_^/9t)] 
n=l 


so 

T'. 


a  E  [a  ,/a 1 1 , 

n=i  '•  ^ 


(Space) 


(Time)  (22) 


Replacing  r,  z  by  r^,  z^  and  equating  these  two  series  term  by  term 

V  =  -2«b£[(a  A„^  ^.j/a  z^)+  (a  A^^  n.i/9t)J  (Space) 

=  -2abjC[a  A^  Jj_j/at] .  (Time)  (23) 

When  this  is  substituted  into  Eq.  (12)  and  the  result  substituted  into  Eq.  (14) 

(Space) 


(24) 


,  -b 

t+£ 

=  -(2iT)"^a  / 

f 

0  ^ 

t-£  ' 

f(x  )  dz  dr  . 

n  n 

n 

=  -(27i)~^  a  f  J 

0 

N 

O 

1 

An  important  difference  in  the  two  cases  may  now  be  noted.  In  the  Time  Case 
only  A^  ^  is  time  dependent.  Equation  (24)  may  be  simplified  by  taking  the  time 
derivative  outside  the  integral. 


-1 

A„  =  (a/at)  (-2  5r)  /  f  ^  a  „  ^  fK>dz„dr^ 

“  "  0  "  z„-£  ^ 


(Time)  (25t) 


L 
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An  iteration  process  will  .^ow  yield 


=  (a>t«)  (-2.)-%"  / 

0  0  0  z^-f  Zq-jC 


A,  e\  )  f(x„).  .  .f(x  )  d  z, .  .  .d  z  d  r,  d  r_.  ..dr  . 
1,0  12  nl  nl2  n 


(Time) 

(26t) 


Usually  the  time  dependence  of  the  initial  field  can  be  separated  Aj^  q 
Aj  0  <r,z)  gj  ^(t)  to  yield 


A  =  (9'^g/at^)  (-277)-"  /  .  .  ./  f  ^  ^ 

0  0  Zg-f 


A',  rt  f(>7,  )•  .  .f(x  )dz,.  .  .dz  dr,.  .  .dr  . 
1,0  1  n  1  n  1  n 


(Time) 

(2Vt) 


If  gj^  Q(t)  =  e“  then 


0  0  z^-i  Zq-C 


a',  n  f(x, ).  .  .f(x^)  dz  .  .  .dz„  dr,.  .  .dr  . 
1,0  1  n  1  n  1  n 


(Time) 

(28t) 


Thus  each  A^  is  cos  t  or  sin  t  times  an  amplitude. 

In  the  Space  Case  such  a  simplification  is  not  possible.  There  are  z  and  t 
dependent  terms  in  the  integral  and  in  the  limits  which  are  not  to  be  differentiated. 
As  will  be  seen  below,  the  successive  A^  differ  by  more  than  simple  phase  shifts. 
However,  the  differential  equation  (Eq.  (9))  has  been  changed  to  an  infinite  number  of 
integral  equations  (Eqs.  (19)  and  (24)).  The  n^*'  equation  corresponds  to  the  n'^^ 
perturbation.  Aq  is  the  initial  field;  A^  is  the  induced  field;  A2  is  the  first  order 
skin  effect;  and  so  on. 


(2Jr)'^  r'^/2  f(Xn) 

,  b  t+£ 


-(277)  ^  Of  /  f  ir^/rr'  O^^^J  f(x,  )dz^dr^ 


(Space) 

(25s) 


b  t+f 


-(277)  f  J  ^  {r^lrr'  [^^^2  +  (9A2  j^/3t)]  f(x2)dZ2dr2 
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A  Fortran  program  for  the  lEiM  7090  has  been  written  to  evaluate 
(2  7r/o!)*^  (n  =  1,  2,  3)  at  r  =  i,  z  -  0.  Slug  length  and  radius  are  variable  inputs  for 
this  program.  Typical  results  are  shown  in  Figure  4.  As  long  as  a  <  2  Jr  the 
approximation 

A  =  Aq  -r  Aj  +  Ag  +  Ag  (Space)  (26s) 

seems  valid. 

3.  CO.NCLI’SIOXS 

The  solution  as  presented  in  Eq.  (25)  allows  one  to  calculate  the  perturbed 
field  around  a  cylindrical  slug  moving  in  a  magnetic  field.  While  Eq.  (25)  is  based 
on  the  assumption  of  the  Aq  given  there,  it  is  easily  extended  to  certain  other 
cases.  Any  other  initial  field  v/hich  has  azimuthal  independence  and  has  no 
azimuthal  component  may  be  accommodated  by  suitably  changing  Aq.  In  practice 
this  would  be  a  solenoid  or  some  approximation  thereof.  If  the  initial  field,  hence 
Aq  is  time  dependent,  a  time  derivative  must  be  included  in  the  equation  for  A^^. 

A^  =  (2J7)-^  a  [^^Aj  q/Sz^)  +  OA^  Q/at)] 

f(x^)  dZj^  dTj^ . 

It  is  possible  now  to  look  back  and  study  the  difference  between  the  Space  Case 
and  the  Time  Case.  The  difference  first  occurs  in  Eq.  (1)  where  a  v  X  B  term  is 
added  to  the  electric  field  in  the  Space  Case.  While  it  appears  that  one  could 
easily  modify  the  electric  field  to  accommodate  this  term,  in  Eq.  (8)  this  term 
gives  rise  to  a  z  derivative  of  the  potential.  In  Eq.  (14)  the  difference  in  these  two 
cases  by  the  inclusion  of  the  variable  t  in  the  limits  of  the  integral  is  the  Space 
Case.  Comparing  the  Space  solution  in  Figure  4  with  the  Time  solution  in  Eq.  (27t) 
shows  that  two  cases,  which  were  apparently  analogous  in  Eq.  (1),  actually  have 
little  or  nothing  in  common.  The  Time  Case  solution  is  the  sum  of  a  number  of 
sinusoidal  waves  of  phase  0,  ±  tt/  2,  and  if.  The  Space  Case  solution  is  the  sum  of 
some  definitely  nonsinusoidal  functions  of  no  simple  phase  relation. 


(Space) 

(27s) 
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Appendix  A 

SinguloritiC'S 


Before  a  computer  can  be  used  to  solve  Eq.  (25),  the  singular  points  must  be 
considered.  Otherwise  the  computer  may  try  to  evaluate  infinity.  In  Eq.  (12),  it 


will  be  noted  that  a  depends  on  (r  x  ) 
n  n 

without  bound  when  r  x 


-1/2 


and  K^.  The  first  of  these  grows 


h 


•  0,  while  K  grows  without  bound  when  x,.-» 

1  /  2 

of  these  is  easily  handled.  Near  r  =  0,  x„  cxrr.  Thus  (rx„)  '  <x  x„. 

^  n  n  n 

of  the  functions  in  powers  of  x  near  x  =0,  shows  that 
^  n  n 


1.  The  first 
An  expansion 


(2 


X  )  K 
n  n 


2  E  oc  X  ‘ 
n  n 


(Al) 


Thus,  at  these  points  a^^  approaches  zero.  Although  this  is  obvious  from  physical 
considerations,  the  computer  must  be  given  mathematical  reasons.  When  x^  =  1, 
a  true  singularity  is  found. 

cc  In  [4  (1  -  Xj^^/^)'^]  .  (A2) 

This  is  the  point  at  which  r  =  r^,  z  =  z^.  In  Eq.  (10)  the  current  cross  section  v.'as 
assumed  infinitesimal.  This  introduces  a  singularity  into  the  current  density.  In 
the  case  of  Aq  one  could  introduce  the  wire  cross  section  but  the  easiest  solution 
is  to  avoid  the  singular  point.  The  slug  will  never  run  through  the  current-carrying 
element  so  the  true  value  of  Aq  at  that  point  is  unimportant.  In  the  other  a^  this 


A2 


singularity  occurs  in  a  function  which  s  to  be  integrated.  Since  is  proportional 
to  the  logarithm,  it  is  to  be  expected  that  the  integral  of  will  remain  finite. 
Physically  this  is  expected  to  be  true  since  infinite  fields  are  not  found  in  the  slug. 
It  appears  to  be  true  mathematically.  W^hile  the  integrals  of  Eq.  (25s)  have  not  been 
solved  analytically,  the  following  integrals  can  be  solved  and  are  all  found  to  be 
finite: 


A 

f  K  dx 
Jq  n  n 

1 

r  K  X  dx 
Vq  n  n  n 


f  K  X  ^  dx  . 
Vq  n  n  n 


In  each  case  the  approximation 

I  S  <*11*  '‘'•n  '  /„  '  <^1  ’  <*n’  «  <*n*)K.l-< 


gives  excellent  results.  With  e  <  10  the  error  is  less  than  1  percent. 
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